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Abstract 


The  aigration  of  slab,  cylindrical  and  spherical  zones  throu^ 
a  block  of  solid  under  the  influence  of  an  electric  field  gradient  has 
been  analyzed  theoretically.  It  was  determined  that  the  zone  could 
migrate  either  up  or  down  the  field  gradient,  E  ,  depending  upon  the 
magnitude  and  sign  of  (i)  the  effective  ionic  mobility  of  the  solvent 
atoms,  U,  (il)  the  Peltier  coefficient  at  the  solid-melt  interface, 
and  (iii)  the  isgposed  temperature  gradient,  G*. 
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THE  MIGSRATIOH  OF  A  LIQUID  ZONE  THROUGH  A  SOLID;  PART  II 


by 

W.  A.  Tiller 

In  a  companion  paper on  this  subject,  a  theoretical  analysis 
of  zone  migration  under  the  influence  of  a  temperature  gradient  as  the 
driving  force  was  treated.  It  was  found  that  the  zone  migration  rate 
depended  upon  the  diffusional  transport  across  the  zone  at  the  atomic 
kinetics  of  both  the  freezing  and  the  melting  processes.  The  stability 
of  zone  shape  was  fotind  to  depend  strongly  on  the  crystallographic 
orientation  of  the  zone,  the  anisotropy  of  the  atomic  kinetics  etnd  the 
orientation  of  the  temperature  gradient  relative  to  the  zone.  In  the 
present  paper,  the  theoretical  analysis  of*  zone  migration  under  the 
Influence  of  an  electric  field  is  presented. 

Theory 

The  starting  geometrical  and  thermal  configurations  for  the 
zones  in  a  block  of  solid  Is  Identical  with  that  given  In  the  first 
paper. The  same  assumptions  will  be  made  here  as  there  and  the  same 
three  zone  shapes  will  be  considered. Only  Isotropic  atomJc  kinetics 
will  be  considered  for  the  quantitative  calculations.  To  best  Illustrate 
the  origin  of  new  effects  associated  with  the  application  of  an  electric 
field,  we  shall  proceed  in  three  steps:  (i)  isothermal,  zero  current 
(G  ■  0,  I  ■  O),  (il)  isotheimal,  non-zero  current  system  (O  -  0, 
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^  ‘K’  y 

.1  j=-  O)  and  (iii)  non-isothermal,  non-zero  current  system  (G  f'  0, 

I  ^  O).  Here  G  and  I  are  the  bulk  temperature  gradient  and  current 
respectively. 


Case  I:  G  =  0,  I  =0 

* 

Consider  an  electric  field  E  applied  to  the  system  along  the 

/ 

Z-axis  such  that  the  potential  increases  with  .Z.  Let  the  effective 
mobility  of  the  solvent  ions  be  Uj  i.e.  if  a  field  of  E  volts /cm  is 
applied  to  the  liquid  alloy,  the  net  rate  of  movement  of  each  solvent 
atom  (whether  ionized  or  not)  is  given  by  U  E  cm/sec.  The  effect  of 
the  electric  field  is  to  force  the  solvent  atoms  in  the  positive  Z 
direction  if  U  is  positive  and  in  the  negative  Z-direction  if  U  is 
negative.  This  additional  force  alters  the  conventional  diffusion 
equation  in  the  zone  so  that  it  takes  the  form 


V 


UE  Sc  1  ^ 
D  SZ  “  D  dt 


(1) 


where  E  is  the  electric  field  gradient  inside  the  zone  and  D  is  the 
diffusion  coefficient  of  the  solvent  atoms  in  the  melt.  Assuming  the 
same  dielectric  constant  outside  the  zone  and  constant  parameters  in 
the  two  phases,  the  electric  field  gradient,  E,  may  be  related  to  the 
bulk  electric  field  gradient,  E  ,  by  the  following  relationships : 


-  3  - 


(i) 

Slab 

* 

E  =  E 

(2) 

(li.) 

Cylinder 

* 

(3) 

(iii) 

Sphere 

3^3  * 

W 

where  ^  la  the  dielectric  constant. 

Since  the  migration  of  the  zone  occurs  at  such  a  slow  rate, 
the  solute  field  in  the  zone  may  he  described,  at  any  moment  of  time, 
by  the  steady-ptate  solution  of  eq.  1,  i.e., 

C  «  A  exp  (UEZ/d)  +  B  (5) 

where  A  and  B  are  constants  (slowly  varying  functions  of  time)  whose 
time-dependence  may  be  determined  by  substitution  in  eq.  1.  For 
positive  U,  the  solvent  concentration  increases  with  Z  and  produces 
zone  movement  down  the  field  gradient..  For  negative  U,  the  solvent 
concentration  decreases  with  Z  and  produces  zone  movement  u£  the  field 
gradient.  In  general,  UEi/D«  1  and  eq.  ^  may  be  used  in  its  linear 
approximat ion . 

The  motion  of  the  zone  destroys  the  coniplete  isothermality 
of  the  system  since  heat  is  evolved  at  the  freezing  interface  and  absorbed 
at  the  melting  Interface.  The  temperature  distribution  in  the  system 
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will  be  similar  to  that  illustrated  in  Fig.  1.  Since  the  tenqperatiire 
distribution  in  the  solid  can  be  shown  to  fall  off  extremely  slowly 
with  distance  from  the  zone  surface,  the  temperature  gradient  created 
across  the  zone  to  provide  the  necessary  heat  transport  is  many  orders 
of  magnitude  greater  than  that  outside  of  the  zone  and  is  almost  exactly 
given  by  setting  G  =  0  in  eqs.  3,  6  and  7  of  reference  1;  i.e., 


(i) 

Slab 

(6) 

(ii) 

Cylinder 

(7) 

(iii) 

Sphere 

(8) 

where  V^,  M  and  K  are  the  velocity  of  the  0  *  TT  interface,  the  latent 
heat  of  fusion  per  unit  volme  and  the  thermal  conductivity  respectively. 
If  the  zone  moves  in  the  positive  Z-direction,  is  positive  and  G 
is  negative  tending  to  Impede  the  zone  motion.  If  the  zone  moves  in 
the  negative  Z-directlon,  should  be  used  in  place  of  and  G  is 
positive  again  impeding  the  zone  motion. 

The  conservation  of  solute  conditions  at  the  zone  Interfaces, 
for  positive  U,  become 


\  K  K  -  “t) 


(9) 
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and 

\  -  “o)  'o  ■  -  ''“““o  -  I  (=<,  -  V’  <“> 

Since  C  -  C  has  the  same  sign  as  UE  and,  since  the  R.H.S.  of  eqs.  9 
O  TT  ^  r 

and  19  have  the  same  and  opposite  sign  as  UE  respectively,  and 

alvays  have  the  same  and  opposite  sign  to  UE  respectively,  i.e.,  if  UE  is 

* 

positive,  is  positive  and  is  negative.  For  UE  positive, 

* 

and  C  <  C  . 

TT  r 

Assuming  that  the  atomic  kinetics  remain  unaltered  hy  the 
presence  and  magnitude  of  E,  eqs.  1,  2  and  3  of  reference  1  may  he  used 
for  the  three  predominsint  growth  mechanisms  considered.  For  the  uniform 
interface  advance  mechanism  we  have,  for  positive  UE,  the  zone  migration 
rate 


[UE/(1  -  k^)] 


.  el  ( +  i) 

can  im  V 

O  TT 


(11) 


>diere  e'  =  -  d/(1  -  k^)  C^,  and  a  =  K^,  Kg  +  and 

2Kg  +  Kj^  for  a  slah,  cylinder  and  sphere  respectively.  For  positive  UE, 
is  positive  and  for  negative  UE,  is  positive.  In  the  denominator 
of  eq.  11,  the  term  in  a  arises  as  a  consequence  of  the  temperature 
gradient  within  the  zone  tending  to  slow  the  motion  of  the  zone.  In  soiife 
systems  ^/a-'lO^  8ec-°C/cm^  and  may  be  as  large  as  the  atomic  kinetic 
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term  except  for  very  small  values  of  Z.  The  appropriate  equations  for 
the  other  mechanisms  may  be  readily  developed  from  the  comparison  of 
eq..  11  with  eqs.  12,  13  and  ik  of  reference  1.  We  can  see  from  eq.  11 
that,  for  positive  E,  the  sign  of  U  deteimines  the  direction  of  zone 
motion. 

Case  II:  G*  =  0,  I*  0 

When  a  current  of  density  I  is  passed  through  the  sample,  a 
quantity,  of  heat  Q  watts/cm^  is  absorbed  at  one  interface  and  generated 
at  the  other  by  the  Peltier  effect.  The  Peltier  coefficient  gP^^  =  O./j 
where  J  is  the  current  density  flowing  through  the  interface  from  solid 
to  liquid  and  of  magnitude  given  by  [  J  |  =  1®  Here,  is  the 

electrical  conductivity  of  the  liquid  and  positive  gP^^  signifies  that 
absorption  of  heat  occurs  when  current  flows  from  solid  to  liquid.  In 
this  case,  a  different  temperature  gradient  is  developed  across  the  zone 
than  that  given  by  eqs.  6,  7  suid  8.  A  term  of  magnitude  +  Q/o£  must  be 
added  to  the  E.H.S.  of  these  equations  (Q  =  •JjE  ^  effects  due 

to  Thomson  heat  generation  in  the  body  of  the  solid  and  the  zone  have 
been  neglected.  Incorporating  the  Peltier  heat  generation  effect,  eq.  11 
beconies 

/  The  value  of  E  is  given  by  eqs.  5,  6  and  7  ln,.reference  1  by  replacing. 

G,  G*  and  K  by  E,  E*  and  a  respectively  with  ^  =  0.  ; 
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L  S  L 

Q!  m 
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Im 


(12) 


When  P  is  of  the  same  sign  as  U,  the  Peltier  heating  tends  to  impede 

the  motion  of  the  zone;  when  they  are  of  opposite  sign,  the  rate  of  zone 

motion  is  enhanced.  Depending  upon  the  magnitude  and  sign  of  P  ,  the 

o  L 

zone  may  move  in  the  opposite  direction  to  that  expected  from  the  sign 
of  UE  alone. 


Case  III;  d*  ^  0,  I*  0 

In  this  case,  we  have  for  the  three  different  atomic  mechanisms, 


-E 

U  e' 

(1-k  )  0!  m 

.  -  0 _ _ _ 

(j  P  +  e£_ 

LS^L  E 

■  e'A  H  €'  p  ' 

am  im  ' 

i.  -  JL 

1  11^ 


(ii)  Screw  dislocation 

(iii)  Two-dimensional  nucleation 


with  (i)  Ifaiform  interface  advance 


(13) 

(IM 

(15) 


r-  V  I 

and  with  p  « 

In  general  e'A  '&/o,  m  <  -  1  and  the  atomic  kinetic  effects  are  again  seen 


,  2K^,  and  3i^  for  slab,  cylinder  and  sphere  respectively. 
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to  depend  upon  i  as  in  the  first  paper. In  figure  2,  V^/s  is  plotted 

versus  u/(l-k^)  for  several  values  of  the  parameter  j6  =  -  (e'/oi  m) 

(Ul  q^-j^  +  PG  /e)  and  we  can  see  the  additive  effect  of  varying  G  and 

different  „P_  on  V^/e.  We  can  see  that  the  direction  of  zone  motion 
b  jj  O 

depends  not  only  on  the  sign  of  E  and  U,  but  also  on  the  magnitude  and 
sign  of  fl). 

Interface  Stability 

We  can  expect  similar  anisotropic  kinetic  effects  as  discussed 
earlier  for  the  zero  electric  field  case.^^^  However,  it  is  likely  that, 
in  some  materials,  the  growth  of  perturbations  of  surface  contour  will 
be  enhanced  while  in  others  it  will  be  retarded  by  the  presence  of  an 
electric  field.  For  a  rumpled  surface  contour,  the  equipotential  lines 
will  no  longer  be  planar  in  the  zone  and  the  current  density  will  also 
be  a  function  of  position.  The  differential  equations  of  heat  and  mass 
transport  and  the  interface  boundary  conditions  will  be  much  more  complex 
than  those  considered  thus  far  due  to  this  positional  variation  of  E  and  J. 
The  solution  to  this  stability  question  will,  in  itself,  be  extremely  com¬ 
plex  and  will  not  be  considered  further  in  this  paper. 

Discussion 

For  some  materials,  as  the  zone  approaches  a  surface  between 
two  phases,  it  will  begin  to  be  influenced  by  the  field  due  to  the  elec¬ 
tric  double  layer  of  this  surface  or  some  other  charge  excess  located 


there.  Such  electrostatic  fields  combine  to  either  enhance  or  retard 

the  motion  of  the  zone  depending  upon  the  strength  and  sign  of  the  field 

* 

and  the  sign  of  the  ionic  mobility.  This  is  a  case  of  J  ■  0  and  -will 

be  more  dominant  the  smaller  are  the  electrical  conductivity  of  the  solid 

outside  the  zone  and  of  the  liquid  in  the  zone.  For  semiconducting  and 

Insulating  solids,  the  width  of  the  normal  space  charge  region  at  a 

“5  -1 

surface  varies  from  about  IG  cm  to  10  cm.  It  seems  possible  that 

such  an  effect  as  this  may  be  responsible  for  the  enhanced  migration 

(2) 

rate  of  Wemick's  '  Ge-Al  wire  zones  through  a  block  of  Ge  as  the  zone 

approached  the  upper  surface.  He  found  that  V  increased  by  a  factor  of 

almost  10  as  the  zone  moved  from  the  lower  to  the  upper  surface  of  the 

block  (temperature  difference  *  30°C).  A  similar  type  of  surface  field 

effect  has  been  observed  for  the  migration  of  spheres  of  aqueous  solution 

(3) 

through  ice  as  they  approached  the  ice-water  interface. 

One  difficulty  with  this  type  of  explanation  for  the  Ge-Al 

wire  is  that  we  normally  expect  the  surface  space  charge  region  in  Ge 

to  be  considerably  smaller  than  the  wire  sizes  used.  Thus,  for  this 

system,  such  an  effect  should  be  observed  only  much  closer  to  the  sur¬ 

fs) 

face  than  observed  by  Wemick'  '  unless  our  estimate  of  the  space  charge 
region  is  small  by  an  order  of  magnitude.  An  altermtive  explanation 
for  this  effect  is  that,  for  cylindrical  and  spherical  zones,  the  value 
of  G  increases  considerably  as  the  zone  becomes  close  to  the  upper  sur¬ 
face  provided  (G  decreases  for  »  K^^).  The  greater  is  the 
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ratio  the  further  from  the  -surface  will  this  effect  of  increasing  G 

begin  to  be  felt.  This  effect  is  illustrated  in  figure  3  for  ^ 

The  isotherms  distort  around  the  zone  to  decrease  the  thermal  gradient 
in  the  zone.  However,  as  this  isotherm  distortion  region  touches  the 
surface,  the  radiation  sees  a  lowered  temperature  and  greater  heat  flux 
enters  this  region  of  surface  to  maintain  it  at  constant  teii5>erature . 

Thus,  the  quantity  of  heat  flowing  into  the  zone  increases  and  G  increases. 
The  magnitude  of  this  effect  can  be  readily  determined  by  considering  the 
electrostatic  analogue. 

The  potential  outside  of  a  conducting  sphere  or  cylinder  Imbedded 
in  a  Tiniform  field,  E^,  is  given  by  that  due  to  a  dipole  or  line  dipole 
respectively  of  moment  proportional  to  E^  and  the  volume  of  the  sphere  or 
cross-sectional  area  of  the  cylinder  respectively.  '  The  field  due  to  a 
pair  of  oppositely  signed  dipoles  aligned  in  the  Z-direction  and  separated 
by  a  distance  2R  produces  a  constant  potential  on  the  median  plane,  distant 
R  from  each.  Likewise,  the  temperature  is  constant  on  the  median  plane 
separating  two  point  doublets  or  line  doublets  of  opposite  sign.  The 
thermal  field  due  to  such  a  pair  of  doublets  immersed  in  a  constant 
temperature  gradient  and  each  situated  distance  R  from  the  external  sur¬ 
face  should  be  the  same  as  that  due  to  a  spherical  or  cylindrical  zone 
placed  distance  R  from  the  surface.  The  temperature  gradient,  G,  in  the 
zone  is  given,  for  these  two  cases,  by; 
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(i)  Sphere 
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(17) 


(il)  Cylinder 


1  + 
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-  1 
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(18) 


We  can  see  that,  as  s/i  - >  «,  eq,s.  17  and  l8  hecome  eq.s.  7  and  6  of 

reference  1  with  V  set  eq,ual  to  zero.  However,  even  for  the  case 

cylinder  and  sphere  respectively.  Even  at  R  only  slightly  greater  than 
i/2,  G  increases  only  505(  and  25^(  for  the  cylinder  and  sphere  respectively. 

Perhaps  the  comhination  of  the  increase  of  G,  increase  of  /, 
faster  atomic  kinetics  and  a  nonzero  value  of  E,  combine  to  give  this 
enhanced  value  of  V. 

When  an  external  electric  field  is  applied  to  the  solvent  block 
such  that  current  J  flows,  the  zone  may  move  either  up  or  down  the  elec¬ 
tric  field  depending  upon  the  parameters  of  the  system  presented  in  eq.  13. 
A  change  in  sign  of  the  ionic  mobility,  U,  can  cause  a  change  in  the  direc¬ 
tion  of  zone  motion  provided  the  first  term  in  the  numerator  of  eq.  13 
dominates.  Similarly,  if  the  second  term  dominates,  a  change  in  sign  of 


K_,  R  ■  /  gives  an  Increase  in  G  of  only  5^  and  2-l/2^  for  the 


P  can  cause  a  change  in  the  direction  of  zone  motion.  One  of  these 
D  L 
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two  effects  is  undoulitedly  responsible  for  the  observations  of  Pfann 

et  al'^'  who  showed  that,  although  J  was  constant  and  was  expected 
'  b  Xj 

to  be  of  the  same  sign  for  a  system  of  binary  alloys  with  Ge,  the  zones 
containing  either  Pt,  Pd,  Au,  or  Ni  moved  in  one  direction  and  the  zones 

containing  either  Ag,  Al,  Bi,  Cu,  Ga,  Sn,  Mg,  Fb,  Sb,  or  Te  moved  in  the 

* 

opposite  direction.  In  this  experiment  G  *  0  so  that  if  U  is  imchanged 
from  alloy  to  alloy,  ®ust  be  of  the  same  sign  as  U  but  with  the  former 
four  alloys  having  a  larger  value  of  -P^  than  the  latter.  Ifowever,  it  is 

b  L 

likely  that  U  will  change  from  alloy  to  alloy  so  that  no  clear-out  conclu¬ 
sion  can  be  made  from  these  observations. 

From  eq..  13,  we  can  see  that  by  varying  G  /e  in  an  experiment  we 

can  determine  the  value.  (6  /e)  that  makes  ■  0  and  thus  can  determine 

'  '0  o 

a  relationship  between  the  constants  in  the  numerator,  l.e.- 


®  ^  -  tr  P 

)c  p 

_€'(l-k^)  L  S^L 

(19) 


From  the  slope  of  the  V^/e  vs  G*/e  curve,  at  large  £  so  that  the  denom¬ 
inator  in  eq..  13  is  ^lnity,  the  value  of  e'p/a  m  may  also  be  determined. 
In  any  electric  field  experiment  with  J  0  we  are  unable  to  separate 

U  and  „P.r  •  Only  if  a  J  ■  0  experiment  can  be.  designed  can  the  term 
o  L 

containing  ^P^  be  set  to  zero  and  a  value  for  U  determined. 

b  Xi 


Since  U  is  the  effective  mobility  for  the  solvent,  it  may  be 
produced  by  the  field- induced  migration  of  solute  ions  since  the  dis¬ 
placement  of  solute  in  one  direction  requires  a  displacement  of  solvent 
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in  the  opposite  direction.  If  a  second  solute  species  is  added  to  the 
melt,  U  will  probably  change  in  magnitude  and  may  even  change  in  sign. 

Conclusions 

The  migration  of  slab,  cylindrical  and  spherical  zones  tinder 
the  influence  of  an  electric  field  has  been  treated.  It  was  determined 

that  the  zone  could  migrate  both  up  or  down  the  applied  field  gradient, 

* 

E  ,  depending  upon  the  magnitude  and  sign  of  the  effective  ionic  mobility 
of  the  solvent  atoms,  U,  the  sign  and  magnitude  of  the  Peltier  coefficient 
at  the  solid-melt  interface,  P  ,  and  the  sign  and  magnitude  of  the  imposed 

D  L 

* 

temperature  gradient,  G  .  It  was  further  determined  that,  in  some  systems, 
electrostatic  effects  could  affect  the  zone  migration  rate.  Finally,  it 
was  shown  that  an  enhanced  migration  rate  of  cylindrical  and  spherical 
zones  in  the  Immediate  vicinity  of  a  surface  will  occur  when  1 

due  to  an  increase  in  G  across  the  zone. 
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FIGURE  CAPTIOUS 


Figure  1 

Figure  2 
Figure  3 


Temperature  distribution  developed  In  ain  Initially  Isothexnal 
system  due  to  zone  migration  under  the  ijtfluence  of  an  Internal 
electric  field. 

Plot  of  Vq/E  veri«i^  U/(l“k  )  for  several  iralues  of  the 
parameter  <fi  •  -  ^  °Pj^  +  pG*/®). 

Isothenn  distortion  due  to  a  sjAierlcal  or  cylindrical  zone  of 
greater  thermal  conductivity  than  the  matrix  vhen  the  zone  Is 
either  imbedded  ixt  the  bulk  or  is  near  a  surface . 
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